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We propose a modified RKKY (Rudermann, Kittel, Kasuya, Yosida) technique to evaluate the 
magnetic properties of the ferromagnetic Kondo-lattice model. Together with a previously developed 
self-energy approach to the conduction-electron part of the model, we obtain a closed system of 
equations which can be solved self-consistently. The results allow us to study the conditions for 
ferromagnetism with respect to the band occupation n, the interband exchange coupling J and 
the temperature T. Ferromagnetism appears for relatively low electron (hole) densities, while it 
is excluded around half-filling (n = 1). For small J the conventional RKKY theory (~ J 2 ) is 
reproduced, but with strong deviations for very moderate exchange couplings. For not-too-small 
n a critical J c is needed to produce ferromagnetism with a finite Curie temperature To, which 
increases with J, then running into a kind of saturation, in order to fall off again and disappear 
above an upper critical exchange J. Spin waves show a uniform softening with rising temperature, 
and a nonuniform behaviour as functions of n and J. The disappearance of ferromagnetism when 
varying T, J and n is uniquely connected to the stiffness constant D becoming negative. 

PACS numbers: 71.10.Fd, 75.30.MB,75.30Vn 



I. INTRODUCTION 
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The Kondo-lattice model (KLMpBHoBOLlEl describes 
the interplay of itinerant electrons in a partially filled 
conduction band with quantum-mechanical spins (mag- 
netic moments) localized at certain lattice sites. Typical 
model properties result from an interband exchange be- 
tween the two subsystems. They are provided by a model 
Hamiltonian, which consists of two parts: 



H = H<+H 



sf- 



(1) 



H s describes uncorrelated electrons in a nondegenerate 
energy band (" s electrons" ) 



n C icr C jcn 



(2) 



where cj (T (ci (T ) creates (annihilates) an electron specified 
by the lower indices. The index i refers to a lattice site 
Hi, and a =f , j is the spin projection. are the hop- 
ping integrals. The second term in Eq. ([!]) represents 
an interband exchange as an intra-atomic interaction be- 
tween the itinerant electron spin cr i and the localized spin 
S t ("/ electrons"): 



H sf = -J 



Si. 



(3) 



According to the sign of the exchange coupling J, a par- 
allel (J > 0) or an antiparallel (J < 0) alignment of itin- 
erant and localized spin is favored with remarkable dif- 
ferences in the physical properties. We restrict our con- 
siderations in the following to the J > case, sometimes 
referred to as " ferromagnetic Kondo-lattice model" , and 
also known as "s — f (s — d) model". The applications 



of the KLM are manifold. In the 1970s it was exten- 
sively used to describe the electronic and magneto-optic 
properties of magnetic semiconductors such as the ferro- 
magnetic EuO and EuSElH Most studies focussed on the 
spectacular temperature dependence of the unoccupied 
band states. The well-known "red shift" of the optic 
absorption edge upon cooling from T = Tq to T = Kl 
can be understood as a corresponding temperature shift 
of the lower bd conduction-band edge. RecenL_quasipar- 
ticle band-structure calculation for bulk EuOEj confirm 
the observed striking temperature dependence as a con- 
sequence of a ferromagnetic exchange J of the order of 
some tenth an eV. The interest in EuO, in particular in 
thin films, was strongly revived by the observation of a 
temperature driven metal-insulator transition below Tq 
with a large drop of— the resistivity by as much as eight 
orders of magnitudcO. This gives rise to a colossal mag- 
netoresistance (CMR) much stronger than that of the 
manganites as Lai-^Sr^MnOa. Recently—a surface half- 
metal-insulator transition was predictedt3 for a Eu(100) 
film. The reason for this is an exchange-split surface 
state. 

In ferromagnetic local moment metals like Gd a 
RKKY-type interaction is believed to cause the ferro- 
magnetic order, so that in principle a self-consistent de- 
scription of magnetic and electronic properties by the 
KLM is possible. The T = 0-moment of Gd is found 
to be 7.63 ^tB"- 7 /zb stems from the exactly half-filled 
4/ shell. The excess moment of 0.63 /iB is due to an 
induced spin polarization of the "a priori" nonmagnetic 
(5d, 6s) conduction-bands, indicating a weak or interme- 
diate J > (H. 

An interesting application of the KLM refers to semi- 
magnetic semiconductors, where randomly distributed 
Mn 2+ or Fe 2+ ions provide localized magnetic moments 
which influence, via the exchange coupling J, the band 
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states of systems like Gai-^Mn^AsEl. The Mn 2+ ions 
provide both localized spins (S = I, concentration x) and 
free carriers (holes, concentration x* < x). The latter 
take care for an indirect coupling between the moments 
leading for small x even to a ferromagnetic order. The 
highest Tc so far reached is 110 K for x = 0.053113 where 
x* is found to be about 15% of the Mn concentration x. 
The origin of ferromagnetism in Gai-^Mn^As is contro- 
versial. The authors of Ref. [l(] claimed that the measured 
Tq{x) is consistent with a conventional RKKY interac- 
tion between the Mn ions. However, they estimated the 
exchange coupling J to be of several eV, therewith com- 
parable to the Fermi energy. This makes second order 
perturbation theory with respect to J rather question- 
able. RKKY should be modified to account for higher- 
order spin polarization terms. In this papar we are going 
to present a respective proposal. 

Another modern application of the J > KLM 
aims at the manganese oxides with perovskite struc- 
tures Ti^D.MnOs (T=La, Pr, Nd; D=Sr,Ca, Ba, Pb) 
which have att racted great scientific interest because of 
their CMREZrEa. Parent compounds like the protagonist 
La 3+ Mn 3+ 03 are antiferromagnetic insulators with a lo- 
calized Mn spin of S = 2. Replacing a trivalent La 3+ ion 
by a divalent earth-alkali ion (Ca 2+ ) leads to a homogen- 
uous valence mixture of the manganese ion (Mnj^, 
Mn 4+ ). The three 3d — t2 g electrons of Mn 4+ can be con- 
sidered as almost localized, forming an S = | spin. Fur- 
thermore, (1 — x) electron per Mn site are of 3d — e g type, 
and itinerant. They are exchange coupled to the S = | 
spins, realizing just the situation which is modelled by the 
KLM. The exchange coupling is ferromagnetic (J > 0). 
Since the manganites are bad electrical conductors, it can 
be assumed that the intraatomic exchange J > is much 
larger than the hopping-matrix .element \t\. With esti- 
mated e g bandwidths of 1 — 2 eVlij and lower limits for J 
of > 1 cVO, the manganites surely belong to the group 
of strongly coupled materials. Many fascinating features 
of the CMR materials can be traced back to an intimate 
correlation between magnetic and electronic components. 
One of the not yet fully understood features concerns 
the spin dynamics in the ferromagnetic phase (x ~ 0.3 
in Lai_ x (Ca,Pb)^Mn03). It was reported 2 !! that the 
spin wave dispersion (SWD) of materials with high Tc 
(e.g.: La .7Pb .3MnO 3 : T c = 355K) can be reasonably 
reproduced by a simple Heisenberg model with nearest- 
neighbour exchange, only. FurukawaEJ gave qualitative 
arguments that this can be understood for strongly ex- 
change coupled materials. However, strong deviations 
from the typical Heisenberg behavior of the .SWD .have 
been found for manganites with lower T c 'sP 2 tP 3 H 24 l It 
is a challenging question whether or not the softening 
of the SWD at the zone boundary can be explained by 
a better treatment of the KLM, i.e. by an appropriate 
consideration of the influence of the conduction-electron 
self-energy on the effecti«e-jCoupling between the loca-L 
ized magnetic momentsE3E3. Solovyev and TerakuraEZI 
argue in such a way claiming that the SWD-softening 



has ajaurely magnetic background. By contrast, Mancini 
et alM believe that the KLM [Eq. (§)] has to be ex- 
tended by a direct (superexchange) interaction between 
the Mn moments to account for antiferromagnetic ten- 
dencies, which dominate the magnetic behaviour of the 
parent material LaMnOa. 

The above-presented, of course incomplete, list docu- 
ments the rich variety of, at least, qualitative applications 
for the KLM. However, the model Hamiltonian [Eq. (jl])] 
provokes such a complicated many-body problem that 
approximations have to be tolerated. Most of the re- 
cent theoretical work on the KLM— aiming at the CMR- 
materials , assumed classical spinsEjLJLJ, mainly in order 
to apply "dynamical mean field theory" (DMFT). It is 
surely not unfair to consider this assumption as rather 
problematic. 

It is the intention of the present paper to use a spe- 
cial RKKY technique for a detailed investigation of the 
magnetic properties of the KLM with ferromagnetic in- 
terband exchange. We exploit the method of Ref. [54] to 
map the interaction operator H s f [Eq. (||)] to the Heisen- 
berg model Hamiltonian. This is done by averaging out 
the conduction electron degrees of freedom by properly 
defined "restricted" Green functions. The resulting ef- 
fective Heisenberg exchange integrals are functionals of 
the itinerant-electron self-energy, and therewith strongly 
temperature and band occupation dependent. These de- 
pendencies manifest themselves in fundamental magnetic 
properties such as phase diagrams, Curie temperature, 
spin- wave dispersions and so on. 

The electronic self-energy is taken from a "moment 
conserving decoupling approach" (MCDA) introduced in 
Ref. |34| and-^viG£6£is£ully applied to several topics in previ- 



ous papers£±]t±lO'LJ. The method carefully interpolates 
between nontrivial limiting cases. 

The final goal of our effort is to obtain quantita- 
tive information about magnetism and the temperature- 
dependent electronic structure of real local-moment sys- 
tems. For this purpose it is intended in a future work to 
combine our model study with a "first principles" band- 
structure calculation,. First attempts in this direction 
were done for EuCo and Geo; however, the magnetic 
ordering could not yet be derived self-consistently. The 
paper is organized as follow. In Sec. |l| we briefly for- 
mulate the many-body problem of the KLM. Section III 
reproduces the most important aspects of the "modi- 
fied RKKY- interaction" (M-RKKY), the understanding 
of which is vital for a correct interpretation of the results, 
which are discussed in Sec. IV. 



II. KONDO-LATTICE MODEL 

The model Hamiltonian [Eq. |l])] provokes a nontrivial 
many-body problem, which can rigorously be solved only 
for a small number of limiting cases. For practical rea- 
sons, it sometimes appears convenient to use the second 
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quantized form of the exchange interaction [Eq. (B)], 

H >f = -\jYS.*rSt nia + S?c\_ a c ia ), (4) 

ia 

with the abbreviations: 

n ia = c\ G c ia , z a = 8 a] - S aL , Sf = Sf + iz a Sf. (5) 

The first term in Eq. (Q) describes an Ising-like interac- 
tion between the z components of the localized and the 
itinerant spin. The second term comprises spin-exchange 
processes between the two subsystems. 

To study the conduction-electron properties we use the 
single-electron Green function, 



we can calculate the spin-dependent occupation numbers 



G ija {E) = ((c irT -J ja )) E . 
Its equation of motion reads 

^ (ESi m — Ti, n ) G m j a (E) — USij — 

I a 

— -^J{z a Iii,j<j{E) + Fu i j rr (E)}, 

where the two types of interaction terms in Eq. (| 
causing the "spinflip function" 



F, 



[E) = {{S7 a c m _ a ;c\ a )) 1 



and the "Ising function" 



(E) = ((Sfc ma ;cl)) 1 



(G) 



(7) 



arc 



(8) 



(9) 



These two "higher" Green functions prevent a direct so- 
lution of the equation of motion (Q). A formal solution of 
the Fourier-transformed single-electron Green function, 



defines the complex electron self-energy Y>\t a (E): 

h 



(10) 



G ka (E) = 



E + ^i + ^o+-e(k)-J: ]ia {E) , 



e(k) are the Bloch energies: 

£ ( k ) = ^E T ^ (R ^ Rj) - 



(ii) 



(12) 



(13) 



-hoc 

r)= j dE f_(E) Ptr (E), 



(15) 



important to fix the band polarization m = (n-f) — 
,f-{E) = (exp(P(E — fi) + l)^ 1 is the Fermi function. 

The central quantity is the self-energy, the determi- 
nation of which solves the problem. However, for finite 
temperatures and arbitrary band occupations an exact 
expression is not available. In Ref. [54| a "moment con- 
serving decoupling approach" (MCDA) predicts the fol- 
lowing structure of the self-energy. 



^ka(E) = -\jz a {S Z ) 



1 



J 2 D*a(E). (16) 



The first term is linear in the coupling J, and propor- 
tional to the the local-moment magnetization (S z ). It 
represents the result of a mean-field approach which is 
correct in the weak coupling limit. The second term 
is predominantly determined by spin-exchange processes 
between band electrons and localized moments. It is a 
complicated functional of the self-energy itself. Thus 
Eq. ([jl]) is not an analytical solution at all, but an 
implicit equation for Y. kcr (E). The MCDA, the de- 
tail of which are presented in ref. Q is a nonpertur- 
bational decoupling approach to a set of properly de- 
fined Green functions that correctly reproduces the non- 
trivial exact limiting cases of the KLM. A weighty ex- 
ample is the special jjase_pf a " ferromagnetically satu- 
rated semiconductor"a'E3'E3. The various " higher" Green 
functions have been approximated by linear combina- 
tion of "lower" functions, where the rigorous spectral 
representations of these "higher" Green functions and 
their exact special cases (S — 1/2, ferromagnetic sat- 
uration, ... ) justify the respective ansatz. Free pa- 
rameters are finally fitted to exactly known spectral mo- 
ments. In a certain sense, the method can be consid- 
ered an interpolation scheme between important limiting 
cases. It WjaS|-SJjte|Cj^sf-|illy used in the past for several 
applicationstjOySESEa. We refer the reader for tech- 
nical details to one of these papers. 

The term D^ a (E) in Eq. ( |l6|) contains certain expec- 
tations values such as 

(a) (n ia ), 

(b) {S?), (SfSf), <(S*) 2 >, 



By use of the quasiparticle density of states 



(c) A CT = {SfnJ, la = (STcU-a), ■ ■ ■ ■ 

(14) By use of the spectral theorem, terms (a) and (c) can 
exactly be expressed by the Green functions (0) , (0) and 
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@. Besides Eq. (]l|) it holds 

1* = -j^ J dE f_(E) lmF iiM {E - M ), (17) 

— OO 
+ 00 

& a = -±- J dE f_(E) ImT ii>il7 (E - /i). (18) 



It remains to express the pure "local-moment" correla- 
tions (b) in terms of the electronic self-energy Ek CT (£0 to 
get a closed system of equations that can be solved self- 
consistently. A finite coupling between the spin operators 
must be of indirect nature since the model Hamiltonian 
[Eq. (|l|)] does not contain any direct exchange. To obtain 
the pure spin correlations (b) , we map the exchange inter- 
action (^) on an effective Heisenberg model with effective 
exchange integrals being functional of the electronic 
self-energy. The result is a "modified" RKKY interaction 
(M-RKKY) that takes into account the exchange- induced 
spin polarization of the band electrons. In lowest order 
it reproduces the conventional RKKY method. 



III. EFFECTIVE EXCHANGE OPERATOR 

The starting idea is to map the interband (s — /) ex- 
change on an effective spin Hamiltonian of the Heisenberg 
type. 



(19) 



For this purpose we use the s — f interaction in the folow- 
ing form: 

H °f = - J 4 E E e_4q Rl ( S * • 4 +qff <W- (20) 



N 

iaa' kq 

a is the band electron spin operator, the components 
of which are Pauli spin matrices. The above-mentioned 
mapping occurs by averaging out the band electron de- 
grees of freedom: 



H 



(21) 



Averaging only in the band electron subspace means that 
(H s f)^ retains operator character in the /-spin sub- 
space. According to Eq. ( p0| ) we have to calculate 



•k+qer C k(T 



q<7 



(22) 



H' has exactly the same structure as the KLM Hamilto- 
nian H [Eq. (pi)] , except for the fact that for the averag- 
ing procedure the /-spin operators are to be considered 
as c-numbers, therefore not affecting the trace. 3' is the 
corresponding grand partition function. The expectation 
value [Eq. does not necessarily vanish for q ^ and 



for a ^ a', as it would do when averaging in the full 
Hilbert space of the KLM. To calculate Eq. ( f22| ) we in- 
troduce a proper "restricted" Green function. 



^k.k+q^) — (( C ka''l c k+q<r))y> 



(23) 



which has the "normal" definition of a retarted Green 
function, only the averages have to be done in the Hilbert 
space of H'. The equation of motion is readily derived 
for both cases, namely, in the case that Cko-' is the active 
operator or t^. + . By use of the "free" Green function, 



G£\E) 



E + n- e(k) ' 
we can combine both equations: 



Gt Z± n (E) — 5 CTCr ' <?q,o (J3)- 



(24) 



r k,k+q 

J 



(25) 



{e-^'>K>G£\E) (S, • -xW, G^(E) 

l( k'- ( k + q)).R lG (o) q (s . . <w' (E) } 



+e~ 



This equation is still exact and can be iterated up to any 
desired accuracy. 

In first order, where the restricted Green function on 
the right-hand side of [Eq. (p5|)] is replaced by the "free" 
Green function Eq. (k4), one obtains 



( G k,k+q(^) 



(i) 



Sera' Sq.Q- 



(26) 



Exploiting the spectral theorem in the restricted Hilbert 
space, for the expectation value [Eq. (p2|)] one finds 

(jj E( c U.^) (c) ) = ^.o- x: /-(e(k)) 

-^E^^-^^U- (27) 

i 

Here we have defined 

D^ = ~ I dEf_{E)lm^Y, G *\ E -ti G< £U E -ti- 

k 



Because of 



a _ 1 / Sf S, 



S : 



(28) 



(29) 



one realizes that the first term in Eq. ( p7| ) does not con- 
tribute to Eq. d2l]), while the second term gives rise to 
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an effective Hamiltonian like Eq. 
effective exchange integrals: 



9J) with the following 



J (1) (q) = 



--J 2 £>W 
2 i 



(30) 



(31) 



E 



/_( £ (k + q))-/_( £ (k)) 
e(k + q) - e(k) 



To obtain the magnetic properties of the local-moment 
system from the effective operator [Eq. (|l9j)], we use the 
spin Green function 



(37) 



is a real number, which 



first proposed by Caller 
eventually helps to derive several spin-correlation func- 
tions by applying the spectral theorem to P^ (E) and a 
proper differentiation with respect to a. The equation of 
motion of p/- (E) reads 



This is the classical RKKY exchange, which comes out 
in first order of our Green function procedure. 

It is well-known that Eq. ( |3l| ) can be equivalently de- 
rived by conventional second-order perturbation theory 
starting from an unpolarized conduction electron gas. In 
order to incorporate the exchange-induced spin polariza- 
tion of the conduction electrons to higher order, it almost 
suggests itself to modify approximation (^) by replac- 
ing in the exact expression [Eq. (p5f)] the restricted Green 
function by the full, all polarization pjocesses containing 
single-electron Green function [Eq. 



GiC\E) 



<5<T"<T^k'k+qGk+qcr(-E), (32) 

8cr>a"Skk>G ktT i(E). (33) 



For the required expectation value [Eq. (|22j)], we obtain 
an expression similar to Eq. (p7|), only Dq has to be 
replaced by 



— OO K 

x (Gf(£-/x) G k+qCT (£- M )- 



(34) 



+4+ q (£-M) gw(£-m) 

Exploiting the inversion symmetry, one proves 

Dn = D$ = ±Y, D «- (35) 

a 

This guarantees an isotropic effective exchange operator 
Eq. (na) with the exchange-integral 



EP\f{E) = 2hS tJ {A a ) - 2?iJ2 J» 



X (((^m^i 

Here we have abbreviated 



(A a ) = ({St,e aS *Sr}-) 



(38) 
(39) 



It is well-known that a simple random-phase- 
approximation decoupling of the higher Green function 
on the right-hand side of Eq. ( |38| ) yields surprisingly 
convincing results in the low- and high- temperature 
region ("Tyablikov-approximation"). Doing so, after a 
Fourier-transformation we obtain 



P^(E) 



2h 2 (A a ) 



E-E(q)+i0^ 



(40) 



In this approximation the magnon energies E(q) are real: 



E(q) = 2n(S z }(Jo-J(q)), 
J = J(q = 0) = ^4 



(41) 



(42) 



Following the Callen-methodEj for the local-moment 
magnetization one obtains 



<S 2 > = h 



(l + S + y) ^ 2S+1 + (S + lf ) (l_+y) 
(1 + tp) 2S + 1 - ip 2S + 1 



2S+1 



ip is the average magnon number, 

^) = ^E e ^ (q 1 )_ 1 - 



(43) 



(44) 



i(q) 



4 



-j-oo 

| dEf_(E)±Yl (~ ^(G^(E-n): 



x G k+q(T (£; - /x) 



(36) 



Via G\ L+<ia {E— n) the effective exchange is a functional of 
the electronic self-energy Eq. (|l6|) , and therefore receives 
a distinct temperature and electron density dependence. 



which determines (S z ) and many other spin correlations. 
Some typical examples are 

(S-S+) = 2ti(S*) <p(S), (45) 

((s*) 3 ) = h 3 s{s + iMS) + h 2 (s z ) (s(s + i)+ 

+^S))-h{(S z ) 2 )(l + MS)). (46) 



Via Eq. ( p6[ ) all spin correlations [Eqs. (|43[)-(|46D] are hxed 
by the electron self-energy [Eq. ([l6|)]. 
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The key quantity of ferromagnetism is the Curie tem- 
perature Tq. Performing the limiting process 



T — ► T, 



(-). 



(s z ) 



H 



one finds, from Eq. m 



k B T c = -h'S(S + l) 



l^(io-J(q) 



(47) 



(48) 



The effective exchange integrals are temperature- 
dependent and have to be used here for T — * Xc- 

For the results presented in Sec. [V, we have used the 
self-energy approach of Ref . j54 . The same approach has 
been applied in Refs. |T^,|l2|ijj3^,||, so that we can re- 
frain here from presenting details of the method. As 
explained after Eq. (|l^) the self-energy of Ref. |34| de- 
pends on three types of expectation values: purely elec- 
tronic terms, which are accessible with the single-particle 
Green function [Eq. ( |q)1 ; pure local-spin correlation such 
as Eqs. (|43|), (f45f), (p6[); and mixed itinerant-electron- 
local-moment correlations, which are expressible by the 
"higher" Green functions [Eqs. (||) and (^)j. Eventually, 
we have a closed system of equations that can be solved 
self-consistently for the desired electronic and magnetic 
properties of the KLM. Some typical results are discussed 
in Sec. Evl. 



IV. MAGNETIC PROPERTIES 

We have evaluated our theory for a s.c. lattice (band- 
width: W = 1 eV) to find out the magnetic proper- 
ties of the ferromagnetic KLM. The latter are, however, 
strongly correlated with the electronic properties, so that 
they can not be understood without referring to the 
itinerant-electron subsystem. Fig. [l] shows the quasipar- 




FIG. 1: Temperature-dependent quasiparticle density of 
states pa-(E) as function of energy (Kondo- lattice model). Pa- 
rameters: W = 1 eV, S = 7/2, J = 0.2 eV, n = 0.2. For these 
parameters the theory yields: Tc = 232 K. The upper half for 
a ="\, the lower half for a =J_. 

tide density of states (QDOS) p a (E) for a band occupa- 
tion n = 0.2 and a moderate exchange coupling J = 0.2 



eV. The self-consistently determined Curie temperature 
turns out to be Tq — 232K. In the low-temperature, 
ferromagnetic phase the QDOS exhibits a remarkable 
temperature-dependence. For T — > the f spectrum 
becomes rather simple (solid line in the upper half of 
Fig. [I]) since the local-moment system is ferromagneti- 
cally saturated. A f electron has therefore no chance to 
exchange its spin with the localized spins. That means 
that only the Ising-like interaction term in Eq. (|J) really 
works, leading to a rigid shift of about —\JS of p-[{E) 
with respect to the "free" Bloch DOS p (E). The down- 
spin QDOS, however, is much more complicated, because 
a I electron can exchange its spin with the saturated 
localized spins. That can, e.g., be done by emitting a 
magnon ?kj(q), where the electron reverses its spin from 
I to t- Consequently, this part of the spectrum ("scatter- 
ing states") exactly coincides with P](E). However, the 
I electron has another possibility to exchange its spin 
with the perfectly aligned localized spins. It can polarize 
its local spin surrounding by repeated magnon emission 
and reabsorption. This can even lead to a bound state, 
which we call the "magnetic polaron". For the special 
case of a single J, electron in an otherwise empty conduc- 
tion band coupled to a ferromagnetically saturated spin 
system the forniajtkm of the magnetic polaron can rigor- 
ously be shownEZlH (see Fig. I in Ref. 32). Formation 
of a polaron costs more energy than magnon emission. 
Polaron states therefore build the upper part of the J. 
spectrum. 

Magnon emission by the j electron is equivalent to 
magnon absorption of a | electron, with one exception: 
magnon absorption can occur only if there are magnons 
in the system. This is not the case at T = in the fer- 
romagnetic saturation. At finite temperature, however, 
there are magnons to be absorbed by j electrons. Conse- 
quently, scattering states appear in the | spectrum too. 
Because of the possible spin exchange both spin spectra 
occupy exactly the same energy region at all T > 0. It 
is interesting to observe that even for such a moderate 
coupling J = 0.2 eV (W — 1 eV) there opens a gap be- 
tween a lower and an upper quasiparticle subband. In the 
lower subband the electron hops mainly over lattice sites, 
where it orients its spin parallel to the local moment, ei- 
ther without or with preceding spin-flip by magnon emis- 
sion (absorption) . The upper subband consists of polaron 
states with finite lifetimes. Note that such elementary 
processes are not recognizable when for mathematical 
simplicity classical spins (S — > oo) are assumedBSr "IP 1 l. 
The QDOS in Fig. [I] clearly demonstrates that the con- 
duction electrons are not at all fully spin polaxizcd as is 
often used in DMFT treatments of the KLMEEil. That 
is true only in the unphysical limit S" ^ oo. 

Let us now concentrate the rest of the discussion on the 
magnetic properties of the KLM. They are determined 
by the respective behaviour of the effective exchange in- 
tegrals Jij. 

Fig. H shows the distance dependence of Jy as it follows 
from the conventional RKKY [Eq. @]. We recognize 
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the well-known oscillating and long-range behavior. Ac- 
cording to Eq. (|3|) two parameters influence the oscilla- 
tion: the interband exchange J and the band occupation 
n. The "amplitude" of the oscillation is proportional to 
J 2 , while the "period" is fixed via the chemical potential 
fi by the conduction electron density n. A remarkable 
temperature dependence is not observed. 

The oscillations of the effective exchange integrals Jij 
are less regular when higher-order polarization effects in 
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the polarized itinerant electron system. Of similar com- 
plexity is the n dependence. In the conventional RKKY 
[Eq. (]3l])] it comes into play only through the chemi- 
cal potential in the Fermi function, while in the full ex- 
pression [Eq. (pq)], several correlation functions as (n a ), 
(S*n a ), (S-°c&_ a ), . . . [Eqs. ©, @, ©] contribute 
to the n dependence. 




R„ (100) 



FIG. 4: The same as in Fig. tt, but for J = 0.2 eV. 



FIG. 2: Conventional RKKY exchange integrals Jij as a func- 
tion of the distance Rij between lattice sites in (100) direc- 
tion for various band occupation n and different interband 
exchange couplings J. Parameters: T — OK, W = 1 eV, s.c. 
lattice. For comparison a respective example of the modified 
RKKY model is included. 



the conduction band are taken into account [Eq. ([56])]. 
Fig. H] and Fig. ^ show two examples for two different 
J's and several band occupations n. The nonregular 
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FIG. 3: Effective exchange integrals Jij of the "modified" 
RKKY model as a function of the distance Rij in (100) di- 



rection for various band occupations. Parameters: 
J = 0.1 eV, W = 1 eV, s.c. lattice, S = 7/2. 



T = OK, 



J dependence results from the fact, that J appears in 
Eq. ([36]) not only as a prefactor J 2 but also in a compli- 
cated manner via the Green function Gk+ qo -(i? — n) of 



To demonstrate the influence of the interband ex- 
change coupling J on the effective Heisenberg-exchangc 
integrals in more detail, in Fig. |] we plot the J de- 
pendence of J x , x = 1, . . . , 9. J x is the effective exchange 
integral between sites separated along the (100) direction 
by x lattice constants. ELI. J4, . . . , Jg show an oscillatory 



0.004 




J[eV] 



FIG. 5: Effective exchange integrals J x between a;th nearest 
neighbors (x — 1,2,... ,9) in (100) direction (s.c. lattice) as 
a function of the interband exchange coupling J. Parameters: 
T = K, n = 0.2, W = 1 eV, S = 7/2. Note the different 
energy scales. 

structure in the weak-coupling region, being, one order of 
magnitude smaller than J2,3 and even two orders of mag- 
nitude smaller than Ji. For J > 0.2 eV J3, . . . , Jg are so 
strongly damped that the magnetism of the KLM can be 
described by a Heisenberg model with next- and next- 
nearest-neighbor exchanges only. In the weak-coupling 
regime, however, the higher terms J3, . . . , Jg can not be 
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neglected. Similar statements can be found vdB. Maybe 
this is an explanation why CMR materials with high 
Tc can be described reasonably welEJ by a short-range 
Heisenberg model, while . oth ep-piaterials with lower Tc 
exhibit strong deviationsE3E§Eil. In any case the main 
contribution to the exchange coupling stems from the 
next-nearest-neighbors. J\ determines the stable mag- 
netic configuration. The negative slope of J\ for J > 0.5 
eV indicates already a breakdown of the ferromagnetic 
order, which is more carefully inspected below. 

The effective exchange integrals J x exhibit an interest- 
ing electron-density dependence (Fig. ^|). For low densi- 
ties up to n — 0.3 only 3\ is of importance. A Heisen- 
berg model with nearest-neighbour exchange will appro- 
priately describe the magnetic properties of the KLM. 
For larger n, however, the exchange interaction becomes 
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FIG. 6: The same effective exchange integrals as in Fig. |5| 
but now as a function of the band occupation n for J = 0.2 
eV. 

very long-range, even Jg is not negligible. From Fig. ^ 
and ^| we learn that the interplay between local interband 
exchange J and band occupation n does lead to a rather 
complicated behaviour of the effective RKKY exchange 
integrals. 

The effective exchange integrals determine the mag- 
netic properties of the KLM. On the other hand, they 
depend on the electron self-energy which is strongly in- 
fluenced by magnetic correlation functions. Typical ex- 
amples of the latter are plotted in Fig. 0, calculated for 
n = 0.2 and J = 0.2 eV. The electronic part of the self- 
consistent solution, represented by the temperature de- 
pendent QDOS, is shown in Fig. The Curie tempera- 
ture is found as Tc = 232 K. The local-moment magneti- 
zation (S z ) is of Brillouin function type, and the qualita- 
tive shapes of the other spin correlations are also familiar 
from the pure Heisenberg model. 

The key-quantity of ferromagnetism is the Curie tem- 
perature Tc- Tc > comes out as a consequence of an 
indirect coupling between the local moments, mediated 
by spin polarization of itinerant electrons. Therefore, a 
strong particle density dependence of Tc has to be ex- 
pected. Figure || shows Tq{u) for various J from the 
weak to moderate coupling regime, calculated according 




T[K] 



FIG. 7: Several spin correlation functions of the local-moment 
system in dependence on the temperature T. Parameters: 
J = 0.2 eV, n = 0.2, S = 7/2, s.c. lattice, W = 1 eV. The 
theory, presented in the text, self-consistently yields Tc = 232 
K. 



to [Eq. (Eg)]. Ferromagnetism appears for low particle 
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FIG. 8: Curie temperature as a function of the band occu- 
pation n for various interband exchange couplings (S = 7/2, 
s.c. lattice, W = 1 eV), calculated by use of the "modified" 
RKKY. model. 

(hole) densities, where the ferromagnetic region increases 
with increasing J. The Tc values are of realistic order. 
A similar n dependence of Tc has been found in Ref. |3^ 
within an extended multi-band KLM. No ferromagnetism 
appears around half- filling (n = 1). It can be speculated 
that antiferromagnetism, which is not considered within 
our approach, becomes stable near n = 1. It is interesting 
to compare the results of Fig. H with their conterparts 
from conventional RKKY (Fig. O), using the same model 
parameters and the same Tc-formula [Eq. (§)]. The Tc 
values are higher, but the region of n where ferromag- 
netism appears, is distinctly narrower. The critical n, 
where Tc vanishes again, is independent of J, in contrast 
to the full theory in Fig. ||. 

The Curie temperature shows a remarkable J depen- 
dence, which is plotted in Fig. |l^ for various band oc- 
cupations. According to conventional RKKY theory one 
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FIG. 9: The same as in Fig. ^, but calculated by use of the 
"conventional" RKKY model. 



would expect a monotonic increase of Tq with increas- 
ing J. A mean-field evaluation of the effective Heisen- 
berg model would lead to To ~ J 2 . Our theory be- 
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FIG. 10: Curie temperature as a function of the interband 
exchange coupling J for various band occupations n (S = 7/2, 
s.c. lattice, W — 1 eV). 



comes identical to conventional RKKY theory for small 
J. However, because of the random-phase-approximation 
(RPA) treatment of the Heisenberg operator [Eq. (fl9|)], 
there are deviations of Tq from the J 2 behavior even for 
small J. A more important and very characteristic fea- 
ture of our modified RKKY approach is the appearence 
of a critical coupling J c for band occupations n > 0.13, 
below which no ferromagnetism occurs. This occurs for 
band fillings, for which conventional RKKY theory fails 
to yield a ferromagnetic RPA solution (Fig. ||). There- 
fore the existence of J c does not disprove the agreement 
of conventional and modified RKKY theory for small J. 
After a steep increase as function of J, Tq runs into a 
kind of saturation, where the plateau becomes smaller 
the larger the bandfilling n. With a further increase of 
J, the Curie temperature again decreases and disappears 
at an upper critical J value. The upper critical value of J 
coincides with the point where the stiffness constant D of 
the spin- wave dispersion [Eq. (|4lJ)] (E(q) ~ Z?q 2 near the 
r point) becomes negative due to a respective behavior 



of the dominating effective exchange integrals J\ , J2 , J3 
(see inset in Fig. || and Eq. (ft8|)). It can be suspected 
that for larger J the local-moment system becomes anti- 
ferromagnetic. Furthermore, it cannot be excluded that 
for still larger J's the system returns to ferromagnetism. 
These features are far beyond conventional RKKY theory 
and due to higher order polarization terms, which come 
into play via the full Green function Gk+qer {E — fi) in the 
expression (|3^) for the effective exchange integral J(q). 
The manganites are .often treated as ferro mag nets with 
J -> 00 in the KLMO. According to Fig. [Hfthis must 
be reexamined. For J rs W ferromagnetism is unlikely 
in the KLM. 
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FIG. 11: Spin-wave dispersion of the ferromagnetic Kondo- 
lattice model as a function of the wave vector for different 
temperatures T. Parameters: J = 0.2 eV, n = 0.2, S = 7/2, 
s.c. lattice, W = leV. The self-consistently calculated Curie 
temperature is Tc = 232 K. 

The spin-wave dispersion, plotted in Fig. |ll| for the 
main symmetry direction, obtains, via the effective 
exchange integrals a distinct temperature dependence. 
Upon heating the dispersion relation uniformly softens, 
disappearing above Tq. This agrees qualitatively with 
the a neutron-scattering study of-,the spin dynamics in 
the manganite Pro.63Sro.37Mn03cHl. However, we do not 
share the view of the authors of Ref. [2^, i.e., that the 
unexcepted results rule out a simple Heisenberg Hamilto- 
nian. In our opinion the Heisenberg model works as long 
as the exchange integrals are renormalized in a proper 
way by the conduction-electron self-energy. 

Changing the doping in the III-V based diluted mag- 
netic semiconductors such as Gai-^Mn^As, or in man- 
ganites like Lai-^Ca^MnOa, does alter the band occu- 
pation and therewith the physics of these materials. It is 
therefore worthwhile to inspect the n dependence of the 
spin-wave dispersion. This is done in Fig. ^ for low ban- 
doccupations and in Fig. |l3] for higher band occupations 
for a model system with J = 0.2 eV. For low densities 
(Fig. |lj) we observe a strengthening of the magnon dis- 
persion with increasing band filling, while the opposite 
is true for stronger fillings (Fig. [l|). This is in accor- 
dance with the Tc behavior of Fig. ||. The band occu- 
pation n = 0.485 is very close to the upper critical point 
of the J = 0.2 eV curve in Fig. ||. The corresponding 
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FIG. 12: The same as in Fig. [H], but for different band occu- 
pations at T = K. 



spin- wave dispersion in Fig. [L5| shows a strong soften- 
ing. A slightly higher n leads to parts of the Brillouin 
zone with negative magnon energies, i.e. the ferromag- 
netic ground state becomes unstable. The stiffness D, de- 
fined by E(q) « D q 2 for q — * 0, becomes negative. The 
points where D = D(n,J,T = 0) = + exactly coincide 
in our theory with those for which Tc = Tc(n, J) — > 0, 
directly derived from the local-moment magnetization 
(S z ). Qualitatively the same n dependence of the spin- 
wave dispersion is reported in Ref. |l]. 



0.05 




FIG. 13: The same as in Fig. |l2| but for higher band occu- 
pations. 

An example for the J dependence of the spin- wave dis- 
persion is presented in Fig.[lJ (n = 0.2, T = K, S = ~). 
Again we observe a strengthening of the dispersion with J 
when Tq increases simultaneously, and a softening, when 
Tc decreases (see Fig. |Io| ). 



V. SUMMARY 

We have evaluated an approximate but self-consistent 
theory for a local-moment ferromagnet described within 
the framework of the ferromagnetic Kondo-lattice model. 



Candidates for this model are magnetic semiconduc- 
tors (EuO, EuS), diluted magnetic semiconductors 
(Gai_ x Mn x As), magnetic metals (Gd, Dy, Tb) and CMR 
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FIG. 14: The same as in Fig. |ll|, but for different interband 
exchange couplings J at T — K. 



materials (Lai^Ca^MnOs). We hapi used a previously 
developed Green function techniques for a detailed in- 
vestigation of the magnetic properties of the exchange- 
coupled local moment-itinerant electron system. An ex- 
tended RKKY mechanism has been worked out to derive 
magnetic phase diagrams as well as spin-wave disper- 
sions. The latter show strong dependencies on tempera- 
ture (uniform softening for T — * Tc), on the band occu- 
pation n, and on the exchange coupling J. These depen- 
dencies are due to the respective behaviors of the effec- 
tive exchange integrals between the localized moments. 
The effective exchange integrals are found by mapping 
the interband exchange (s — /) coupling, characteristic 
of the KLM, to an effective Heisenberg model. They 
turn out to be functions of the electronic self-energy, be- 
ing therewith strongly temperature and electron density 
dependent. For small J our theory reproduces the results 
of the well-known conventional RKKY treatment, which 
represents second-order perturbation theory to the KLM. 
Already for very moderate J, however, substantial devi- 
ations appear, when higher-order terms of the induced 
conduction-electron spin polarization bring their influ- 
ence to bear. 

For the future the presented approach has to be ex- 
tended to antiferromagnetic moment configurations in 
order to complete the phase diagram. Furthermore, 
Coulomb correlations have to be introduced for the 
conduction-electron subsystem. Their neglect in the orig- 
inal KLM seems to be only poorly justified. 

It is further intended to apply the modified RKKY pro- 
cedure to the antiferromagnetic KLM (J < 0) in order to 
investigate the interplay between " Kondo screening" and 
RKKY. This, however, requires a seripus check whether 
or not the MCDA in its present forrnEJ is able to account 
for the subtle low-temperature (Kondo) physics. 
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